Error— proof programmable self— assembly of DNA— nanoparticle clusters 



Nicholas A. Licata, Alexei V. Tkachenko 
Department of Physics and Michigan Center for Theoretical Physics, 
University of Michigan, 450 Church Str., Ann Arbor, Michigan 48109 

We study theoretically a new generic scheme of programmable self-assembly of nanoparticles 
into clusters of desired geometry. The problem is motivated by the feasibility of highly selective 
DNA-mediated interactions between colloidal particles. By analyzing both a simple generic model 
and a more realistic description of a DNA-colloidal system, we demonstrate that it is possible to 
suppress the glassy behavior of the system, and to make the self-assembly nearly error-proof. This 
regime requires a combination of stretchable interparticle linkers (e.g. sufficiently long DNA), and 
a soft repulsive potential. The jamming phase diagram and the error probability are computed for 
several types of clusters. The prospects for the experimental implementation of our scheme are also 
discussed. 
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I. INTRODUCTION 

Over recent years, significant attention has been at- 
tracted to the possibility of nanotechnological applica- 
tions of DNA H , , [| , (3 , @] . Among the various pro- 
posals, one of the most interesting directions is the use 
of DNA-mediated interactions for programmable self- 
assembly of nanoparticle structures [1] , Q j Hi > • Sev- 
eral schemes of such self-assembly have been studied both 
experimentally and theoretically. Their common theme 
is the use of colloidal particles functionalized with spe- 
cially designed ssDNA (markers), whose sequence de- 
fines the particle type. In such systems, selective type- 
dependent interactions can be introduced either by mak- 
ing the markers complementary to each other, or by us- 
ing linker-DNA chains whose ends are complementary to 
particular marker sequences. 

Recent theoretical studies have addressed the expected 
phase behavior [pcL melting properties [TTj ]. and mor- 
phological diversity [l2| of DNA-colloidal assemblies. In 
particular, there are indications that these techniques can 
be utilized for fabrication of photonic band gap materi- 
als [H[j[3- Despite significant experimental progress, 
the long-term potential of DNA-based self-assembly is 
far from being realized. For instance, most of the exper- 
imental studies of DNA-colloidal systems report random 
aggregation of the particles [lj|. Some degree of struc- 
tural control in these systems has been achieved, mostly 
by varying the relative sizes of particles, rather than by 
tuning the interactions 

In the present paper, we take a broader view of pro- 
grammable self-assembly. While this theoretical study 
is strongly motivated by the prospects of DNA-colloidal 
systems, our main objective is to address a more general 
question: how well can a desired structure be encoded 
by tunable interactions between its constituents? The 
particular model system on which we focus consists of 
distinguishable particles with individually controlled in- 
teractions between any pair of them. In the first stage 
of our study, we analyze a simplified yet generic version 
of such a system. All the particles have the same re- 



pulsive potential, while the attraction is introduced only 
between selected pairs of particles in the form of a spring- 
like quadratic potential. This potential mimics the effect 
of a stretchable DNA molecule whose ends can selectively 
adsorb to the particular pair of particles. We then in- 
troduce a number of additional features which make the 
model a more realistic description of an actual DNA- 
colloidal system. 

Our major result is that a combination of stretchable 
interparticle linkers (e.g. sufficiently long DNA), and a 
soft repulsive potential greatly reduces (or totally elim- 
inates) the probability of self-assembling an undesired 
structure. The experimental prototype is a system of 
particles in a mixture of two types of DNA molecules 
which can selectively adsorb to the particle surface. The 
first type are DNA molecules with two "sticky" ends, i.e. 
both end sequences of the DNA are complementary to 
the particle marker sequence. With one end adsorbed 
to the particle surface, the remaining sticky end makes 
it possible to introduce an attractive interparticle poten- 
tial between selected particle pairs. The second type are 
DNA molecules with one sticky end which adsorbs to the 
particle surface. These DNA strands give rise to a soft 
repulsive potential of entropic origin between all particle 
pairs. 

There is a natural analogy between our problem and 
the folding of proteins where interactions between amino 
acids encode the overall structure. However, it should be 
emphasized that the task of programmable self-assembly 
in colloidal systems is even more demanding than pro- 
tein folding: any intermediate metastable configuration 
has a much longer lifetime and therefore means a mis- 
folding event. Because of this, we were looking for a 
self-assembly scenario which does not require thermally 
activated escape from a metastable configuration. This 
makes the problem additionally interesting from the the- 
oretical perspective of "jamming" , a phenomenon ac- 
tively studied in the context of granular and colloidal 
systems. Our results can be interpreted in terms of a 
jamming-unjamming transition controlled by the inter- 
action parameters. 

It should be emphasized that the goal of this work 
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is primarily conceptual, as opposed to providing a man- 
ual for the immediate experimental realization of ordered 
colloidal structures. Nevertheless, future experimental 
schemes will be forced to overcome obstacles presented by 
colloidal jamming. With this in mind, one of the most 
salient features of our model is the ability to smooth the 
energy landscape by tuning the interactions between par- 
ticles. 

The plan for the paper is as follows. In section II, we 
address the problem within a simplified generic model 
which mimics the nanoparticle system with stretchable 
DNA connections. An unexpected and very encouraging 
result of this study is that the misfolding (or jamming) in 
the model system can be completely avoided for a certain 
set of parameters. In section III the original model is 
adapted to a more realistic situation which incorporates 
the random character of the DNA-mediated interactions. 
In section IV we discuss the prospects for the future ex- 
perimental implementation of our scheme. In section V 
we summarize the major results. 

II. BEADS AND SPRINGS MODEL 

Consider an isolated group of repulsive particles linked 
via a polymer spring to their desired nearest neighbours. 
We assume that the DNA marker sequence is the same for 
any two markers attached to the same particle, but dif- 
ferent particles have different marker sequences (i.e. each 
particle has a unique code). In this case, the attraction 
between any two particles can be effectively switched on 
by adding DNA "linkers" whose ends are complemen- 
tary to the corresponding marker sequences of the parti- 
cles. As a first approach to the problem, we introduce 
a generic "Beads & Springs" model which incorporates 
essential features of the DNA-nanoparticle system. The 
model system contains N particles with pairwise (type- 
independent) repulsive potential U(r). In general, this 
repulsion may have a hard-core or soft-core behavior, or 
be a combination of the two. In order to model the DNA- 
induced type-dependent attraction, we introduce a har- 
monic potential (linear sp ring s) which acts only between 
selected pairs of particles[18|]- Thus, the model Hamil- 
tonian has the following form: 

H = 2 ^2 kJ <*13 \ r a ~ 173I 2 + U(\r a - 173I) (1) 

a,f3 

Here a, (3 are the particle indices, r a are their current po- 
sitions, and K is the spring constant. The connectivity 
matrix element J a p may be either 1 or 0, depending on 
whether the two particles are connected by a spring, or 
not. Our goal is to program the desired spatial config- 
uration by choosing an appropriate connectivity matrix 
J a p. A natural construction is to put J a p = 1 for any 
pair of particles which must be nearest neighbors in the 
desired cluster, and not to connect the particles otherwise 
(i.e. put J a j3 = 0). This construction assures that the 
target configuration is the ground state of the system. 



Note that our problem is somewhat similar to that of 
heteropolymer folding. In that case, the selective inter- 
actions between monomers (e.g. amino acids in protein) 
are responsible for the coding of the spacial structure of 
a globule. Our major concern is whether the kinetics of 
the system will allow it to reach the ground state within 
a reasonable time. Unfortunately, since the Brownian 
motion of a typical nanoparticle is relatively slow (com- 
pared to molecular time scales), it is unrealistic to expect 
that our system will be able to find the target configura- 
tion by "hopping" between various metastable states, as 
in the case of protein folding. However, our case is dif- 
ferent because the attractive force grows with distance, 
as opposed to the short-range nature of heteropolymer 
self-interactions. As we shall see below, this difference is 
essential, making it possible for the system to reach the 
ground state without stopping at any metastable config- 
uration. 

We have performed a molecular dynamics simulation of 
the above model by numerically integrating its Langevin 
equation: 

b- 1 r a = -V a H + T la (2) 

Here b is the particle mobility. The thermal noise has 
been artificially suppressed in this study (i.e. r) = ). 
In other words, we have assumed the worst case scenario: 
once the system is trapped in a local energy minimum, 
it stays there indefinitely. The equations of motion were 
solved numerically by a first order Runge-Kutta method. 
First, we studied a system of N = 49 distinguishable 
particles in 2D, whose native configuration was a 7 x 7 
square cluster (see Figure[T]). Their initial positions were 
random, and the connectivity matrix was constructed ac- 
cording to the above nearest-neighbor rule. 

First, we studied the case when U(r) is a hard-core 
potential. More precisely, the repulsive force was given 
by a semi-linear form: 

hc(r) = -^^- = Ko(d-r)Q(d-r) (3) 

Here O is the unit step function, and d is the diameter 
of the hard sphere. The parameter kq determines the 
strength of the hard-core repulsion, and it does not affect 
the results, as long as kq k. In our simulations, we 
found that the hard sphere system eventually stops in a 
configuration definitely different from the desired one, a 
behavior which is well known in the context of granular 
and colloidal systems as "jamming" [l9j. Remarkably, 
the jamming can be avoided when the hard-core repulsion 
is replaced by a soft-core potential: 

U sc (r) = U eiq>(-r/\). (4) 

Here the decay length A is of the order of the equilibrium 
interparticle distance r . This indicates that the energy 
landscape can be made smooth by a combination of long- 
range selective attraction and soft-core repulsion. 
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a) Hard-core repulsion b)Soft-shell repulsion 




FIG. 1: (Color Online). Programmable self-assembly in 2D, 
studied within the Beads & Springs model. 49 particles are 
all distinct and connected with linear springs to encode the 
desired configuration (7x7 square). The jamming behavior is 
observed for the case of hard spheres (a). However, the as- 
sembly of the target structure can be achieved if the repulsion 
is sufficiently "soft" (b). 



The result is surprising and remarkably robust. In 
particular, in order to expand our finding to the 3D case, 
we studied the self-assembly of particles into tetrahedra 
of various sizes (N — 10,20,35). This time, the hard 
core interaction potential was superimposed with a soft 
shell repulsion, which makes the model more relevant for 
an actual DNA-colloidal system: 



U (r) = U hc (r) + U sc (r) . (5) 



After the system has fully relaxed, a geometric measure 
of the folding success is determined by comparing particle 
separations of the desired final state to those generated 
from a set of random initial conditions. Figure [3] shows 
the "jamming phase diagram" for these systems. To as- 
sign a point on the diagram to the correct folding regime, 
we required 100 consecutive successful folds. While this 
criterion can only give an upper bound on the jamming 
probability (which is approximately 1%), an additional 
analysis gives strong evidence that the correct folding re- 
gion of the diagram corresponds to zero probability of 
jamming. 
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FIG. 2: (Color Online). "Jamming phase diagram" obtained 
for programmable self-assembly of tetrahedral clusters, within 
the Beads & Springs model. The control parameters depend 
on the equilibrium interparticle distance r , the diameter of 
the hard sphere d, and the range of the soft-shell repulsion A. 



III. SELF-ASSEMBLY IN DNA COLLOIDAL 
SYSTEMS 

As we have seen, the introduction of a soft-core repul- 
sive potential U sc is crucial to a successful self-assembly 
proposal. In a real system, this repulsion can be gener- 
ated e.g. by DNA or another water soluble polymer ad- 
sorbed to the particle surface. The mechanism is quite 
independent of the monomer chemistry, but for the sake 
of concreteness we will speak of the repulsion generated 
by DNA. Namely, we assume that a certain fraction of 
the DNA "arms" of the " octopus- like" particles are not 
terminated by a sticky end, and only play the role of a 
repulsive "buffer". When the polymer coverage is suf- 
ficiently low, the interparticle repulsion is primarily due 
to entropy loss of a chain squeezed between two parti- 
cles. The characteristic length scale of this interaction 
is given by the radius of gyration of the "buffer" chain, 
R g . The corresponding repulsive force f sc can be calcu- 
lated exactly in the limit of relatively short buffer chains, 
R g <C d. The result of this calculation [12| can be ade- 
quately expressed in the following compact form: 

hc() d(r-d) P { 2RI ) W 

Here is the total number of buffer chains per parti- 
cle. Note that even though this result is only valid for 
Rg <C d, it correctly captures the Gaussian decay of the 
repulsive force expected for longer chains as well. There- 
fore, we expect the results to be at least qualitatively 
correct beyond the regime of short buffer chains. 

In addition to the modified soft potential, we have 
to take into account the random character of the real- 
istic DNA-mediated attraction. It originates from the 
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fact that (1) the number of DNA "arms" of the orig- 
inal octopus-like particles will typically be determined 
by a random adsorption process, and (2) the fraction of 
the DNA chains recruited for linking a particular pair 
of particles is also random. In terms of our original 
model, this means that the "springs" will not have the 
same spring constant. If the individual linkers are mod- 
elled by Gaussian chains [2(| , the overall spring constant 
for a particular pair of connected particles is given by 
K a p — kBTm a fj/2Rg , where R g is the radius of gy- 
ration of a single linker, and m a p is the number of in- 
dividual chains connecting the particles. We assume 
that this number obeys the generic Poisson distribution: 
P(m) = m m e~ m /ml. As formulated, the model is cast 
as a system of coupled differential equations: 



--^-^Japm a f3r a f3 + f hc [r a p) + f sc (r a p) 



1q/3 



P 

(7) 

Here r Q(3 = |r a - rg|, n Q/3 = (r Q - 173) /r af3 . 

We have studied the behavior of the system as a func- 
tion of two dimensionless parameters, one of which is the 
ratio of the buffer radius of gyration to the particle di- 
ameter, Rg/d . The other parameter characterizes the 
relative strength of the attractive and repulsive forces: 



m / R, 



(8) 



Including the above modifications of the model, the es- 
sential result of the study is that the jamming probabil- 
ity can be drastically suppressed, similarly to the previ- 
ous case. However, the jamming in this system cannot 
be eliminated completely. Instead of an actual "phase 
boundary" , we have observed a sharp crossover to the 
regime of predominantly good folding, in which the error 
probability is suppressed to a modest level ~ 10 — 20%. 
Interestingly, the behavior is nearly independent of the 
energy parameter a, which makes R g /d the only major 
control parameter. Figure shows the error probability 
as a function of this geometric parameter for tetrahedral 
clusters of different sizes. As this plot indicates, the mis- 
folding behavior gets suppressed as R g /d exceeds 1. 

As a further test of the robustness of the model, we 
consider a modified attractive potential which deviates 
from the linear Hooke's law. For larger forces we enter 
the Pincus regime [2l|, where the end-to-end extension 
of the polymer chain r is related to the external tension 
/ ~ n. This tension law incorporates the excluded 
volume interaction of individual linker DNA with them- 
selves, which was not previously considered. Remark- 
ably, the major result for error suppression carries over, 
as illustrated in figure [U 

Our description has several limitations which are to 
be addressed in future work. In particular, our discus- 
sion is only applicable to the limit of modest coverage of 
particles with buffer chains. This case of weakly overlap- 
ping adsorbed chains (known as the " mushroom regime" ) 
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FIG. 3: (Color Online). Error probability P as a function of 
aspect ratio R g /d for tetrahedral clusters with modified soft- 
potential and realistic DNA-mediated attraction. Each data 
point on the misfolding profile represents 100 trials. 
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FIG. 4: (Color Online). Error probability P as a function 
of aspect ratio R g /d for tetrahedral clusters in the Pincus 
regime. Each data point on the misfolding profile represents 
100 trials. 



is drastically different from the high-coverage "polymer 
brush" behavior [22|. Nevertheless, our major conclu- 
sions appear to be rather robust. The condition for error 
suppression (i? g /(f > 1) is the same in both the harmonic 



and Pincus regimes. 



IV. DISCUSSION 

Previous studies have demonstrated an experimen- 
tal implementation of self-assembly in a DNA-colloidal 
system [23(, 24]. The approach in these studies differs 
from our vision of controllable self-assembly in two ma- 
jor ways: (1) because the "linker" DNA chains(^ 20nm) 
used are much shorter than the particle diameter(~ 
l/zm), the particles behave as sticky spheres, and (2) 
there is little diversity in particle type, where the struc- 
tures result from interactions of a one or two component 
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system. We would like to discuss issues related to a fu- 
ture experimental implementation similar to the modified 
Beads & Springs model. 

In our molecular dynamics simulations we assumed 
that (1) the desired group of particles has already been 
localized in a small region of space, and that (2) the in- 
terparticle connections have already been made. There 
are a number of experimental challenges associated with 
implementing the self-assembly proposal of our simula- 
tions. In another manuscript [25| we provide a detailed 
discussion of the localization problem, which is the first 
major experimental intermediary. 

After localization, the next step is to make the desired 
connections between particles within the cluster. To do 
so one can add short ssDNA with sequences sasb to link 
particles A and B. The DNA marker sequence sa for 
particle A is a sequence of nucleotides complementary 
to the sa portion of the linker sequence sasb- The hy- 
drogen bonding of complementary nucleotides forms base 
pairs which join both marker strands to the linker, cre- 
ating a DNA bridge between the two particles. After 
the interparticle links are formed, they should be made 
permanent by ligation. Since the spring constants of 
the above dsDNA chains are too small to drive the self- 
assembly of a desired cluster, we propose to melt them 
either by changing the temperature or pH. As a result, 
the dsDNA links will be turned into ssDNA with a much 
higher effective spring constant (due to the shorter per- 
sistence length). This will trigger the self-assembly sce- 
nario similar to the one discussed within the Beads & 
Springs model. Note that DNA entanglements may be 
effectively eliminated if the procedure is done in the pres- 
ence of DNA Topoisomerases. 

V. CONCLUSION 

In conclusion, we presented a model of DNA- 
colloidal self-assembly which exhibits a tunable jamming- 



unjamming transition. The combination of a soft-core 
repulsion with a type-dependent long range attraction 
provides a natural tunneling of the energy landscape to 
the ground state configuration. This is to be contrasted 
with the case of protein folding, where under physio- 
logical conditions the interactions between amino acids 
are screened to several angstroms. Because this length- 
scale is much shorter than the spatial extent of the native 
structure, large regions of the energy landscape are flat, 
which prohibits formation of the native state on the basis 
of tunneling alone. As a result the folding rate is nec- 
essarily limited by the diffusion of amino acid seg ments 
looking for their desired nearest neighbors [2q|.[27l| The 
fact that the potential in the DNA-colloidal system is 
long ranged is essential, allowing us to avoid the pitfall 
of slow particle diffusion. 

Within the Beads and Springs Model, we obtained the 
jamming phase diagram for several modest sized tetrahe- 
dral clusters. We identified a regime of parameter space 
with error-proof folding, and demonstrated the impor- 
tance of introducing a soft-core repulsion. The origi- 
nal model was then adapted to include several features 
of realistic DNA-mediated interactions. Although the 
jamming cannot be completely eliminated in the modi- 
fied system, we identified a regime of predominantly good 
folding, and calculated the error probability for tetrahe- 
dral clusters. The jamming behavior is determined by 
a single geometric parameter R g /d. We concluded by 
discussing prospects for an experimental implementation 
of our self-assembly scheme. 
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